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We calculate the inclusive production of h c (ht) at e + e~ colliders with the center-of-mass energy 
from the CLEO-c energy to Z° boson mass at leading order of nonrelativistic QCD. At Z° boson 
mass, the cross sections are 39 ~ 703fb for h c (lp), 37 ~ 61fb for ht(lp) and 44 ~ 73fb for hb(2p). 
At the B-factory, it is 86 ~ 212fb for h c (lp). For h c at the CLEO-c and h b (lp,2p) at the B- 
factory, the perturbative QCD expansion is not good and the results are much smaller than the 
experimental measurements. It is clearly shown in all the results that the color-octet state ( 1 S'o) 
contributes dominantly while the color-singlet state (Pi) contributes small or even negative part, in 
contrast to the case of inclusive J ftp production where the color-singlet state is found to contribute 
dominantly. 

PACS numbers: 12.38.-t, 12.39.St, 13.60.Hb, 14.40.Pq 



In recent years, many experimental measurements for 
P-wave quarkonia h c ,hb(l + ~) have been achieved. The 
related branch ratios were measured [l~6], the masses of 
them are measured precisely 04H|) and the cross sec- 
tions for h c (hb) production via e + e~ annihilation at the 
CLEO-c (B-factory) are also measured [l(| E| ■ More ex- 
perimental measurements on P-wave quarkonium states 
h c ,hb could be expected in the future. It provides a 
new place to test or improve our knowledge of quantum 
chromodynamics (QCD) on heavy quarkonium produc- 
tion and decay. 

To study heavy quarkonium decay and production pro- 
cesses, nonrelativistic QCD(NRQCD) 12] is a successful 
factorization theorem, in which calculation is factorized 
into process-dependent short-distance coefficients, to be 
calculated perturbatively in the strong-coupling constant 
a s expansions, and universal long-distance matrix ele- 
ments (LDMEs), to be extracted from experiment. In 
this scheme, the QQ pair can be produced in any Fock 
state n — 2S + 1 l}^ with definite spin S, orbital angular 
momentum L, total angular momentum J, and a = 1, 8, 
where a = 1 is for color-singlet (CS) and a — 8 is for 
color-octet (CO) states which finally evolve into physical 
quarkonia through nonperturbative processes. The rela- 
tive importance of all the CO and CS states are estimated 
by velocity scaling rules, which weigh each of the LDMEs 
by a definite power of the heavy-quark velocity v in the 
limit d«1. In this way, the theoretical predictions are 
organized in double expansions in a s and v. 

Based on heavy quark spin symmetry of NRQCD, the 
LDMEs for h c (hb) are simply related to that for Xc(Xb), 
therefore, studies in the past already supply the infor- 
mation of the LDMEs for h c (hb) from the fit of Xc(Xb) 
related processes. Experimental data [IH, [lij shows that 
J/tp production from Xc feeddown count on a large part 
in J /jp hadroproduction. Based on the experimental 
data |l3[ and LO calculation on Xc hadroproduction, an 
estimate of the CO LDME for Xc and h c was given in 
Ref. [TH, [l6j]. Employing these matrix elements, the cal- 
culations of h c hadroproduction at the Tevatron [l7| and 



LHC dHH predicted a significant yield. Photoproduc- 
tion of h c was investigated in Ref. [2(| by using a CO 
LDME extracted from the decay B — > XcJ + X, the re- 
sults indicated a significant cross section at the DESY 
HERA. Recently, the Xc hadroproduction calculated up 
to QCD next-to-leading order (NLO) with an estimate of 
the CO LDME for Xc is given in Ref. [Hj]. In the study 
of polarization for prompt J/ip hadroproduction 22j , the 
CO LDME for Xc is given by fitting the Tevatron and 
LHCb data. 

In this paper, we study heavy quarkonium P-wave 
states h c (hb) inclusive production via e + e~ annihilation 
by calculating the related cross sections at the leading- 
order (LO) of NRQCD. Detailed results are given at the 
Z° boson peak, B-factory and CLEO-c. The discussions 
on the comparison between our calculation and the ex- 
perimental measurements are given. 

In NRQCD framework, at LO of a s and v 2 , cross sec- 
tions for H state production can be expressed as 

^W = / lp m(O ff ( 1 P 1 [11 )) + / ls [ 8] (O ff ( 1 4 81 )), (1) 

where f n denotes the short-distance coefficient cor- 
responding to the NRQCD operator O h (n) and the 
LDMEs of the operators are 

(0 H CS [ * ] )) = (0\^T a xa f H a H ^T a m, (2) 

J. J. k 

W|0>, 



(o h Cp[ 1] )) = (o|V t (-* T )x< • wh t : 



where H represents the hadron state h c , or 
cc( l Pf-\ 1 sf ] ) states. The CS operator O ff ( 1 P 1 [1] ) 
is of order v 2 while the CO operator O h ^ sj^) is of 
order v°. However, to hadronize into h c , the CO state 
have to emit at least one soft gluon, which rises the 
order of the LDME (O h " ( 1 S 1 ^)) to a s v 2 , the same order 
of v as the CS LDME (O^pf 1 )). 

The leading processes for inclusive h c production are 
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listed as 



e+e- ^ccCPl 1] )+g + g, 
e+e- -> cc( l P^ ] ) + c + c, 



(3) 



It is easy to see that the three processes are of the same 
order of v 2 at LO. The cross sections for the second and 
third processes are infrared (IR) divergence free and can 
be calculated directly. The cross section for the first pro- 
cess, although at LO, is IR divergent when one of the 
gluons is soft. To deal with the divergence, the NRQCD 
factorization formulas involved are as follows, 

°{h c ) = f lp m (O h «(P[ 1] )) + / lg w (O^ 1 ^ 1 )) (4) 

and 

°es l * ] ) = f 1P w(o ls ° ] epi l] )) +f ls m(o is[ ° ] esW))(5) 

The matrix elements of operators are calculated in di- 
mensional regularization scheme at LO with an ultravio- 
let cutoff (NRQCD scale) /l«a = m c and are given as: 



(0 1 4 81 (ipW)} = o, 



(6) 



where N c is 3 for SU(3) gauge field and u c is defined as 
1 



— -lE + UlnC-^) 
em a Ma 



(7) 



with /i/j being the renormalization scale. The divergence 
in this LDME will cancel that in a (ipf 1 ), which can 
be isolated by using the two-cutoff phase space slicing 
method [231 as 



(8) 



where <r s and <r s are from the gluon-soft and hard region, 
respectively. The boundary of the two regions is E g — 

^5 si where E g is the energy of the soft gluon, \fs is the 
center-of-mass (CM) energy of e + e~ colliders and 6 S is 
an arbitrary positive number small enough to make the 
soft approximation good enough. The soft part can be 
expressed as 



Snmi Ni: ° n 



(9) 



where 
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Using Eq.© and Eq.©, we obtain the expressions of 
short-distance coefficients in Eq.(j4]) as 



A, 



.[8] 
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(11) 

(12) 



where 



s + 4m^ 



e e 6 s — 4m 2 4/i« ^ 

It is clearly shown that all the short-distance coeffi- 
cients are IR divergence free and the cross section for the 
first process is well defined. To calculate ct( 1 S'q 8 ') and 
a" s ( 1 p] 1 ') for the first process and the cross sections for 
the second and third processes, we apply our Feynman 
Diagram Calculation package(FDC) [24| to generate all 
the needed FORTRAN source. 

From the heavy quark spin symmetry of the NRQCD 
Lagrangian, it is obvious that the LDME (O' 1 " ^S^)) 
for the intermediate state cc( 1 S'g 8 ') evolving into h c is ex- 
actly the same as that for the intermediate state cc( 3 S^) 
evolving into \ c at LO of v 2 . It gives 



(o h -Cs l o ] )} « (o Xci ( 3 si 8| )>- 



(14) 



We will employ the LDME value obtained in Ref. [l6[ as 

(0* C1 ( 3 S[ 81 )) = 0.0098 GeV 3 , which is from a fit of the 
Tevatron experimental data with the LO calculation for 
transverse momentum p t distribution of Xcj hadropro- 
duction and is independent on the NRQCD scale n\ 
since the range of pt in the fit does not include p t — 
where /i\ is involved in the calculation to factorize IR 
divergence. However, the result in our calculation is 
//A dependent, therefore we have to fix a value of /xa- 
Since /itA-dependent term is proportional to In(fiA) and 
it is found that different choices of /ia in a reasonable 
range around the heavy quark mass tuq do not cause 
larger uncertainty in final results when the CM energy 
is far above 2toq, the default value ^a = wq is cho- 
sen for h c production at both B-factories and Z-factory, 
and for hb production at the Z-factory. For h c produc- 
tion at the CLEO-c and hb production at B-factories, 
we give results for different fi\ choices to show the un- 
certainties caused by this parameter. In the numerical 
calculation, we have the following common choices as 

(0*"(if ] )) = §K C (0)| 2 = 0.322 GeV 5 , a = 1/128, 
sin 2 6W = 0.226, m z = 91.2 GeV and T z = 2.49 GeV. 
And other parameters are explicitly given in each indi- 
vidual case. The soft cutoff 8 S independence is checked 
in the calculation and S s — 0.0001 is used. 

For the CLEO-c experiment at the CM energy y/s = 
4.17GeV, only the first and last processes in Eq.Q con- 
tribute and the gluons in the first process are too soft 
that the convergence of perturbative calculation in QCD 
is not good. Therefore, it is not surprising to see that 
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the theoretical results can not describe the experimen- 
tal measurement. To show the uncertainty from the 
c-quark mass m c , the QCD coupling constant a s and 
NRQCD scale /za, our results are listed in Table. U by 
choosing m c as 1.3 GeV, 1.5 GeV and 1.76 GeV, a s as 
a s (3.0 GeV) = 0.26 and a s (*/s/2) = 0.30, and fi A as m c 
and m c /2. 



a s (2m c ) - 0.26 


m c ( GeV) 




°o 9 


Total 


1.3 


-4.449(-1.534) 


4.406 


-0.04(2.87) 


1.5 


-3.057(-1.560) 


3.014 


-0.04(1.45) 


1.76 


-1.656(-1.104) 


1.531 


-0.12(0.43) 


otsWs/2) - 0.30 


m c ( GeV) 




°o 9 


Total 


1.3 


-5.923(-2.042) 


5.084 


-0.84(3.04) 


1.5 


-4.071(-2.077) 


3.478 


-0.59(1.40) 


1.76 


-2.204(-1.470) 


1.766 


-0.44(0.30) 



TABLE I: Total Cross sections(pb) for 

l<j[8] 



production from 
a s and 

Values in the brackets correspond to 



Pi 1 ', ' S^' and the sum of them with different m. 
Ha at y/s = 4.170 GeV. 
MA = m c /2. 




FIG. 1: total cross section for production of h c via e eT 
annihilation as a function of CM energy. The parameters are 
chosen as m c = 1.5 GeV and a s = 0.26. 



The results in the table are all negative for (j-A = m c 
and become positive for fj-A = m c /2 while the experimen- 
tal measurements at the CLEO-c [l(| give the cross sec- 
tion for e+e" -> h c (lP)+Tr+n- as 15.6±2.3±1.9±3.0pb. 
On the results, there are three points to be addressed. 

First, we focus on the a s expansion. The value of a s 
depends on the renormalization scale fj,R, and in the total 
cross section, there are logarithm terms ln(E 3 / (j>r) which 
must be small to make sure that they do not ruin the ex- 
pansion. E s can be m c , s and p 9 max = (s — Am 2 )/(2y/s). 
The smallest scale p 9 max i s 1-0(0.6) GeV when m c is 
1.5 GeV(mft, c /2), as a result, there is not a proper choice 
of fj,R to lead to good convergence in perturbative QCD 
expansion. The perturbative expansion becomes better 
when the CM energy is larger than y/s = 6 GeV where 
"Pmax i s 2 GeV. For inclusive h b production at the B- 
factory, it becomes better when the CM energy is larger 
than 12 GeV where p 9 max is 1.8 GeV. 

Second, we investigate the v 2 expansion in shot- 
distance coefficients. The worst propagators for the v 2 
expansion in the first process in Eq.([3|) is 



where pn, p g and q denote the momentum of h c , gluon 
and the relative momentum of c-quarks in h c , respec- 
tively. It is easy to obtain a, a Lorentz invariant, as 
a = \q\cos9q tPg /m c in h c rest frame where \q\/m c is the 
velocity v of c-quark in the quarkonium. Therefore the 
a 2 expansion in the propagator is the nonrelativistic ex- 
pansion. It is easy to see that validity of the expansion 
only depends on the value of a 2 and has no relation with 
the momentum of the gluons. In our case of h c , v 2 is 
thought to be larger than that of J/tp, so the expansion 
is better for J/ip than that for h c . 



Third, we have to check more detail of NRQCD fac- 
torization scheme. The default NRQCD scale fi\ — 
m c sets up a boundary between perturbative and non- 
perturbative parts. The maximum energy of the soft 
gluons in the non-perturbative part is /ia = m c which 
even exceeds the maximum gluon energy p g max = 1 GeV 
in the phase space in the e + e~ — > h c (lP)+gg at CLEO-c 
energy. So /i\ = m c is not a suitable choice, and /iA = 
m c /2, which is smaller than p 9 maxl is chosen to present 
our results. Finally it is clear that the uncertainty from 
the choice of NRQCD scale ha is larger near production 
threshold and with choice of /ia — rn c /2, results for h c in- 
clusive production at the CLEO-c is 0.30 ~ 3.04pb in our 
calculation which is just 1/50 ~ 1/5 of the measurement 
of the exclusive process e + e~ —> h c (lP) + it + it~ . 

In Fig. [TJ we present the dependence of the total cross 
section for the inclusive h c production on the CM en- 
ergy with m c = 1.5 GeV and a s = 0.26. When the CM 
energy is smaller than 6 GeV, the perturbative QCD ex- 
pansion is not good and the total cross section turns out 
to be sensitive to the choice of NRQCD scale fiA and is 
of large uncertainty. The uncertainty in our result be- 
comes smaller quickly as the CM energy goes larger. We 
find that there are two peaks in the plot, one is at Z° 
boson mass and the other is at about 5 GeV. And the 
ratio of the cross sections at the two peaks is about 1.1:2, 
which means that they are in the same order of magni- 
tude. It can also be seen that the e + e~ — > h c (lP) + gg 

and e + e~ — > /i c ( 1 S^) + g dominate the inclusive produc- 
tion rate when the CM energy is smaller than 20 GeV, 
while e + e~ — > h c (lP) + cc dominate that when the CM 
energy is larger than 20 GeV. 

All the above discussions are applicable to hb produc- 
tion at the B-factory with better convergence for non- 
relativistic expansion since v 2 becomes smaller. In the 
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numerical calculation, we employ the LDMEs, of which 
the CO ones [l5[ are estimated based on the NRQCD 
scaling rule and the CS ones |25| are obtained based 
on potential model, as (O h "(lP){p[ 1 ')) = 6.1 GeV 5 , 
{O h »(2P)(p[ 1] )) = 7.1 GeV 5 , (CM^Sf)) = 

0.43 GeV 3 and (O h » (2P)( 1 S l 8] )) = 0.52 GeV 3 . In Ta- 
ble. HQ and HH the results for h b (lP) and h b (2P) pro- 
duction at the B-factory are listed, where [iA = m b /& is 
chosen with the same consideration as for h c case. We 
can see that the uncertainty from the choice of NRQCD 
scale fi\ becomes smaller than that in h c production at 
the CLEO-c. In comparison with the experimental mea- 
surement in Ref. [11] . where exclusive production rate 
of e+e" -> h b n+TT- are 0.416pb and 0.695pb for h b {lP) 
and h b (2P), respectively, our results are only 1/13 ^1/5 
and 1/19 ~ 1/7 for h b (lP) and h b (2P), respectively. 



a a (2m b ) - 0.18 


m b ( GeV) 




s o 9 


Total 


4.75 


-3.77(-1.50) 


76.7 


73(75) 


4.95 


-2.78(-1.34) 


53.1 


50(52) 


5.13 


-1.88(-1.05) 


32.7 


31(32) 


a s (v^/2) -0.21 


m b ( GeV) 


L P[ il 99 


s o 9 


Total 


4.75 


-5.13(-2.04) 


89.5 


84(87) 


4.95 


-3.78(-1.82) 


62.0 


58(60) 


5.13 


-2.56(-1.43) 


38.2 


36(37) 



TABLE II: Total Cross sections(fb) for h b (lP) production 
from 1 Pl 1] , 1 S 1 q ] and the sum of them with different m b , a s 
and ha, at ^/s = 10.865 GeV. Values in the brackets corre- 
spond to fiA = mb/4. 



a 3 (2m b ) - 0.18 


m b { GeV) 


l P?99 


9 


Total 


4.75 


-4.39(-1.75) 


92.8 


88(91) 


4.95 


-3.25(-1.56) 


64.2 


61(63) 


5.13 


-2.19(-1.22) 


39.5 


37(38) 


a s {y/i/2) -0.21 


m b ( GeV) 




s o 9 


Total 


4.75 


-5.98(-2.38) 


108.3 


102(106) 


4.95 


-4.42(-2.12) 


74.9 


70(73) 


5.13 


-2.98(-1.67) 


46.1 


43(44) 



TABLE III: Total Cross sections(fb) for h b (2P) production 
from 1 P1 1] , 1 S l o ] and the sum of them with different m b , a a 
and ha, at \fs = 10.865 GeV. Values in the brackets corre- 
spond to fiA = m b /4:. 



Fig. [2] shows the dependence of the total cross sec- 
tion for inclusive h b (lP) production on the CM energy. 

1 f8l 

We can see that the S process dominates throughout 
the \/~s range, and the first peak where the cross section 
reaches its maximum value is at the CM energy 12 GeV, 
while Z° boson mass is the second peak. And the ratio 
of the cross sections at the two peaks is about 4:3, which 
means they are almost of the same value. 

We calculate the production of h c at the B-factory by 
changing quark mass and QCD coupling constant to show 
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FIG. 2: total cross section for production of h b via e + e _ 
annihilation as a function of CM energy. The parameters are 
chosen as m b — 4.75 GeV and a s = 0.18. 

the uncertainties caused by different choices of the two 
parameters. The results are listed in Table. IIVI In this 
case, theoretical predictions should be better than those 
at the CLEO-c since the uncertainty from the choice of 
//a becomes much smaller. 



a s (2m c ) - 0.26 


m c ( GeV) 


L H%9 




9 


Total 


1.3 


8.4 


46.2 


162.3 


212 


1.5 


-6.5 


15.5 


137.7 


147 


1.76 


-12.2 


3.7 


113.5 


105 


OsWs/2) -0.21 


m c ( GeV) 


l H%9 




°o 9 


Total 


1.3 


5.5 


30.2 


131.1 


167 


1.5 


-4.2 


10.1 


111.2 


117 


1.76 


-7.9 


2.4 


91.7 


86 



TABLE IV: Total Cross sections(fb) for h c production from 
^f 1 , X 5^ 81 and the sum of them with different m c and a s , at 
= 10.6 GeV. 



The results at Z° boson mass are listed in Table. IV; 
Table. Elland Table. EED for h c (lP), h b (lP) and h b (2P) 
respectively. 

Fig. 131 shows the distribution of invariant mass Mx 
of the final states with h c excluded for h c production in 
the CS processes at B-factory and Z-factory CM energy. 
The parameter choice is m c — 1.76 GeV and a s = 0.26. 
Fig-IUgives that for h b (lP) and h b (2P) production. The 
CM energy is 91.2 GeV and b-quark mass for h b (lP) and 
h b {2P) are 4.95 GeV and 5.13 GeV, respectively. QCD 

coupling constant is a s — 0.18. The Sq process con- 
tribute in the phase space point where Mx is zero. How- 

1 rgi 

ever, S Q state evolves into a P-wave quarkonium re- 
quires it emitting soft gluons which will combine with 
the other gluon and hadronize. Therefore, its contribu- 
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a s (2m c ) - 0.26 


772 c (, LreV J 


-fi 99 


/ | cc 


1<j 8 „ 

6 o 


lotal 


1 Q 

1.6 


z.Ol 


fine 
Dyo 


z.Ul 


(Uo 


1.0 


l.oa 


ooy 


1.(4 


o4z 


1. 1 D 


u.oy 


101 


1/18 
1.4o 


lOo 


a s (Vi/2) - 0.13 


m c ( GeV) 


L P^99 


\pWcc 


O o y 


Total 


1.3 


0.63 


175 


1.01 


177 


1.5 


0.38 


85 


0.87 


86 


1.76 


0.22 


38 


0.74 


39 



TABLE V: Total Cross sections(fb) for h c production from 

" !?[8] 



lp[l] 

s/s = 91.2 GeV. 



S [ °> and the sum of them with different m c and a s , at 



a s (2m b ) - 0.18 


m b { GeV) 


^00 


L P[%1 


D o 9 


Total 


4.75 


0.73 


10.3 


49.6 


61 


4.95 


0.63 


8.24 


47.5 


56 


5.13 


0.54 


6.79 


45.8 


53 


a»(V£/2) -0.13 


m b ( GeV) 


l Pl%9 


L P[%b 


'S^9 


Total 


4.75 


0.38 


5.37 


35.8 


42 


4.95 


0.33 


4.30 


34.3 


39 


5.13 


0.28 


3.54 


33.1 


37 



TABLE VI: Total Cross sections(fb) for h b (lP) production 
from 1 P[ 1] , 1 S l o ] and the sum of them with different m b and 
a s , at ^ = 91.2 GeV. 



tion to the Mx distribution is around Mx = 0, but not 
exactly at Mx — 0. In the same situation, the first term 
in Eq. (fT2|) contribute to the same phase space point of 
the two body final states h c + g while the other gluon 
is soft and has already been integrated out. It is in the 
same way to hadronize as the CO process. It is obvi- 
ous that the Mx distribution is not a good observable 
at small Mx in theoretical calculation. Although there 
is no way to do this hadronization and give right distri- 
bution since it is a nonperturbative part, there is a try 
on this problem for J/ip production at e + e~ colliders in 
Ref HI. 

In summary, we study the inclusive production of 
h c (lP), hb(lP) and hb(2P) via e + e~ annihilation with 
the center-of-mass energy from the CLEO-c energy to 



Tota 
h G gg 




M»(GeV) 

FIG. 3: distribution of invariant mass of the final states 
excluding h c for the process e + e~ — > h c X at the B- 
factory (upper) and Z- factory (lower). 




M»(GeV) 

FIG. 4: distribution of invariant mass of the final states ex- 
cluding h b for the process e + e~ — > h b X for h b (lP) (upper) 
and h b (2P) (lower) at the Z-factory. 



a s (2m b ) - 0.18 


m b ( GeV) 


l P[%9 


L P^bb 


l S^9 


Total 


4.75 


0.85 


12.0 


60.0 


73 


4.95 


0.73 


9.61 


57.4 


68 


5.13 


0.63 


7.92 


55.4 


64 


a s (Vi/2) -0.13 


m b ( GeV) 


l P[%9 


1 P[%1 


l S^9 


Total 


4.75 


0.44 


6.26 


43.3 


50 


4.95 


0.38 


5.02 


41.5 


47 


5.13 


0.33 


4.13 


40.0 


44 



TABLE VII: Total Cross sections(fb) for h b (2P) production 
from 1 P[ 1] , 1 S l o ] and the sum of them with different m b and 
a„, at ^ = 91.2 GeV. 



Z° boson mass at leading order of nonrelativistic QCD. 
A detailed discussion on the validity of perturbative cal- 
culation for the total cross section of h c (hb) production at 
the CLEO-c(B-factory) is given. It is easy to see that the 
perturbative QCD expansion is not good in these cases. 
It is also found that validity of the nonrelativistic expan- 
sion only depends on the value of v 2 and has no relation 
with the momentum of the related gluons, so the nonrel- 
ativistic expansion is available. Besides, /ia = m c is not 
a suitable choice, and /i\ — m c /2, which is smaller than 
the maximum gluon energy 1 GeV in the phase space, 
is chosen to present the results. In our numerical calcu- 
lation, the LDMEs (O^fS^)) and (O^^S 1 ^)) from a 
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fit of the Tevatron experimental data with the LO cal- 
culation for transverse momentum p t distribution of Xcj 
hadroproduction are independent on the NRQCD scale 
/i\ , since the range of pt in the fit does not include p t = 
where p\ is involved in the calculation to factorize IR 
divergence. Therefore, in our results, ha independence 
can not be achieved and it can bring large uncertainty 
when the CM energy approach the production threshold. 
At Z° boson peak, the cross sections are 39 ~ 703fb for 
h c (lp), 37 - 61fb for h b (lp) and 44 - 73fb for h b (2p). At 
the B-factory, they are 86 - 212fb for h c (lp), 32 - 87fb 
for hb(lp) and 38 ~ 106fb for hb(2p) where the produc- 
tion rate for hb(lp, 2p) are about 5~19 times smaller than 
the experimental measurements. And at the CLEO-c, it 
is 0.3 ~ 3.1pb, which is of large uncertainty from the 
choice of p\, far away from the experimental measure- 
ment. In another aspect, the results clearly show that 
the color-octet state ( 1 iSq) contributes dominantly while 
the color-singlet ( 1 P 1 1 ) state contributes a small part, 



in contrast to the case of inclusive J/il> production at 
the B-factory, where the color-singlet state is found to 
contribute dominantly. From previous works at NLO of 
QCD for inclusive J/ip production at e + e~ colliders, it 
is clear that the QCD NLO correction can enhance the 
production rate 1.7 times for e + e~ — > J/ip + cc [IS] and 
1.2 times for e + e~ — > J/tp + qq[28{. It is quite interest- 
ing to see the effect from QCD NLO correction on the 
inclusive h c (hb) production in future work. 

Note added When we were preparin g ou r draft, a re- 
lated work appeared at arXiv recently [29| . which study 
on h c production at B-factories. Once take their renor- 
malization scheme in the calculation of NRQCD correc- 
tions to CO operator and their parameters, our calcula- 
tion on h c production agrees with their numerical results. 
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